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Abstract 



> 

, An action functional for the loop quantum cosmology difference equa- 

1 tion is presented. It is shown that by guessing the general form of the 

solution and optimizing the action functional with respect to the param- 
eters in the guessed solution one can obtain approximate solutions which 
are reasonably good. 



o 

CT 1 - 1 Introduction 

u 

One of the main candidates for quantum gravity is loop quantum gravity which 
is canonical quantization of gravity using connection variables. It is background 
independent and non-perturbative. This last property makes it suitable for con- 
sidering quantum effects on black holes and big bang. Applying the theory to 
■ cosmology leads to loop quantum cosmology (LQC)pQ whose main equation for 

the evolution of the scale factor, that is the Hamiltonian constraint, is a differ- 
ence equation. This difference equation may be solved numerically [3J, or exactly 
in some cases [2]. Also continuous approximation, in which the difference equa- 
tion is approximated by a differential equation, maybe useful at least for some 
parts of the domain of the difference equation[IJI|]|. The most important result 
of LQC is the resolution of big bang singularity Also there are phenomeno- 
logical results from the effective equations that incorporates the modifications 
of matter hamiltonian due to quantum gravity ( see e.g. 

Here we shall study the possibility of adopting approximate methods based 
on variational methods for solving LQC difference equation. First we shall 
briefly review the variational methods for difference equations[7j. Then, we 
shall write an action functional appropriate for the difference equation of LQC 
and obtain the first integrals of the equation. Next we shall find approximations 
to the solutions, not by solving the difference equation but by guessing a solution 
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and optimizing the action, and compare these approximate solutions with the 
exact ones. 



2 Review of variational methods for difference 
equations 

Here we summarize the results of ref. [Jj for a case suitable for our problem. 
Suppose that we have an action of the form: 

M 2 

A= C ( n l Qn, Qn-1, Qn-2] Q* n , Ql-v QI-2)) (i) 

n=Mi 

in which Q n is the degree of freedom of the system with Lagrangian C. The 
equation of motion is derived using the least action principle stating that the 
action should be stationary with respect to the variations Q n — > Q n + SQ n . The 
variation in the Lagrangian is: 

5£(6Q n ,6Q* n ) = K (1) SQ n + K {2) 8Q* n (2) 

where the operators TZ^ and TZ^ are defined as: 
2 

=Y, d Qn-^Qn,Qn-l,Qn-2;Qn,Q* n -l,Q* n -2)(L {1) ) k (3) 
fc=0 

2 

=Y / d Q: _ k ^n;Q n ,Q n ^,Q n ^;Ql,Q; i _ 1 ,Q; i _ 2 )(L^) k (4) 

fc=0 

and the lowering operators and are defined as: 

L {1) X(Q n ) = X(Q n _i); and L {2) X(Q* n ) = X(Q* n _ 1 ) (5) 
Defining the adjoint operators: 



11 = X] 5 Q"-2+fc £ ("+ 2 - fc ; Qn+2-k, Qn+l-k, Qn-k; Q* n +2-k, Qn+l-k: Qn-^i^T 
k=0 

(6) 

2 

^ =51 ^Q^_2+fc^(" +2 ~ fc ' Q'l+S-fc) Qn+l-k, Qn-k', Qn+2-ki Qn+l-k, Qn-k)(£ )* 
k=0 

and using the Lagrange identity proved in ref. [7|: 

K {1) SQ n = 5QjZ {1) l+AB^\SQn,Q n ); and Tl {2) 5Q* n = 6Q* n K (2 ' ] l+AB^(5Q n ,Q* n ) 

(8) 
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where A means evaluation of the difference of the quantity at n — 1 and n and 
with: 

2 i 

B« (SQ n , 8Q* n ) = -J2Yl Wn-i^-i 

i=l j=l 

£(n — j + i; Q n -j+i, Qn-j+i-li Qn-j+i-2, Qn-j+i, Qn-j+i-1, Qn-j+i-2) (9) 

2 i 

B™{5Q n ,5Q* n ) = -^C/gv 
i=i j=i 

C(n - j + i; Qn-j+n Qn-j+i-1, Qn-j+i-2, Qn-j+i, Q*i-j+i-l, Qn-j+i-2) (10) 

and substituting in the action, one sees that the A terms are zero by not varying 
the boundary (initial) conditions and the 1Z adjoint terms leads to the equation 
of motion as: 

2 

^Q n C(n + 2 - fc; Qn+2-k, Qn+l-k, Qn-k, Q*n+2-k, Qn+l-k, Qn-k) = i 11 ) 

fc=0 

and its complex conjugate. 

In order to find first integrals of the above difference equation, one first 
should define a symmetry of the Lagrangian. A Lagrangian is said to be sym- 
metric under the transformation: 

Qn^Qn + eu^(n;Q n ,Q n ); Q* n - Q* + eu^(n;Q n ,Q* n ) (12) 

provided one can find some v(n; Q n , Q n -i,Qn-2, Q n , Q n -ii Qn-2) f° r which: 

SC = eAv(n; Q n , Q n _ u Q n - 2 ,Q* n , Qn-nQn-2) (13) 

If a Lagrangian is symmetric in this way, we have a constant of motion as: 

2 

J = ev(n; Q n , Q„_ 1; Q n - 2 ,Q* n , Q* n -i,Q* n -i) - E B (t) (u^(n; Q n , Q* n )) (14) 

3 Variational principle and first integrals for LQC 

The suitable Lagrangian for LQC difference equation is: 

C = K'^Qn - Q n _ 2 | 2 + (k-^Q - 2) + S n _i) |Q n _!| 2 (15) 
leading to the well-known equations of motion: 

Qn-2 — QQn + Qn+2 = KS n Q n (16) 

where S n is the matter contribution to the difference equation. 1 

1 Thc quantum dynamics of LQC states is obtained from the Hamiltonian constraint dif- 
ference equation: 

(v e+5 i - v l+Uo )e- 2lVl ° *i>(4>,e + u ) - n(v e+to - v e - to )f (</>,£)+ 
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Let us now investigate the first integrals. First consider the case f2 = 2 and 
S n = 0. Then the Lagrangian has translational symmetry generated by = 1. 
This leads to v = 0. The corresponding conserved quantity is: 



J — —(Qn+l + Q n - Qn-1 — Qn-2 + Qn+1 + Qn ~ Qn-1 ~ Qn-2) i^V 

which can be simply checked that it is conserved. For this case the equation of 
motion can be solved as Q n = (A + Bn)e %p7rn — constant where p is an integer, 
and A and B are constants [3]. It can be checked that J = constant. 

A more useful result can be obtained when one considers the complete 
Lagrangian (with matter and curvature contributions). The Lagrangian has 
global phase transformation invariance. That is to say, under transformations 
Qn — * Qn& la or in its infinitesimal form Q n — > Q n + iuQ ni the corresponding v 
is zero, and the conserved quantity can be found as: 

J = ia(Q n+1 Q n -i - Q n +iQ* n -i + QnQn-2 - QnQ* n - 2 ) ( 18 ) 

This conserved quantity can give some information about the relation between 
the real part and imaginary parts of Q. Assuming Q n — p n + iq n one shall 
get X n + X n -\ = Ci where C\ is a constant, and X n = p n +iq n -i - Pn-iQn+i- 
The solution for X n is X n = — % + (— l) n_1 C2 where C2 is another constant. 
Writing q n = p n Y n we have Y n+2 —Y n = 21 - ' 11 which can be solved for 
Y n leading to the following relation between the real and imaginary parts of Q n : 

q - = c + c 4 (-ir + \ Ea + (-i)"- i ) c,2( ~ 1)J ~ c ' l/2 (19) 

Pn 2 ^ PjPj+2 

where C3 and C4 are two other constants. 



4 Solving the difference equation 

The variational principle can also be used for deriving approximate solutions 
for cases in which one has not an exact solution. This can be done by guessing 
a solution and choosing the parameters such that the action takes its extremum 
value. We shall illustrate this method with some examples for which we have 
the exact solutions. 

(Vi-M - Vi-wo ) e2ire ° t - m = ~K~f 3 4i 2 p c^(e, 4>M4>, 

in which Q = 2 — 4£g7 2 r(r — 1), C m is the matter (</>) Hamiltonian, and £ p is the Planck 
length, r is the spin connection parameter denned as V = k/2, so f2 = 2 for a flat model and 
SI = 2+^"q7 2 > 2 for a closed model. £o and fi are quantum ambiguities, 7 is Immirizi— Barbero 

parameter, re = 8irG, and Vg = (^^) ^ ^\ are the eigenvalues of the volume operator. 
Changing to the variable Q{t) = (V t+ t - Vt-i )e~ ilr / 2 e- iG "t'ip{<l>, t) and choosing £ = 2£ n 
leads to the difference equation mentioned in the text. Although this discussion is about a 
massless scalar field, it is not hard to see that other matter sources as well as cosmological 
constant leads to the same form of difference equation. 
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5 Case: Q = 2 and S n = 

The action is: 

oo 

A = K- X \Qn~Qn-2? (20) 

71 — — OO 

One can guess a solution of type e an . Substituting in the action we have: 

oo 

A = 4/c" 1 sinh 2 a e {a+a " )n (21) 

n— — oo 

Minimizing this leads to either a = or a = ipn. Both are exact solutions. A 
more general guess for Q is n^e Qn , leading to action: 

oo 

A = k~ 1 J2 e {a+a " )n \e a {n + lfe- a {n-lf\ 2 (22) 

n— — oo 

Minimizing this leads to either a — ipn, (3 = for which the solution is Q n = 
constant, or a — ipn, (3 — 1 where the solution becomes Q n — n. Both are 
exact solutions 



6 Case: Q ^ 2 and S n = 

Now the action would be: 

oo 

A=K~ 1 J2 (lQ«-On-2| 2 + (^-2)|Q„_ 1 | 2 ) (23) 
n— — oo 

Guessing a solution of the form e^ n , we have: 

oo 

first term = n- 1 ^ 4 sinh 2 /?e (,3+ ' 3 '> (24) 

TL— — OO 
OO 

second term = k- 1 (SI - 2)e ( ^+' 3 * )n (25) 

71— — OO 

The problem is that these summations are divergent, so first we have to make 
them meaningful by changing the summation limits to (— N ■ ■ ■ N). Let's assume 
(3 is real and positive. Then we have 

g-2iV/3 _ 20{N+l) 

^l = K- 1 {4sinh 2 /3 + fi-2}- (26) 

1 — e " 

Since A*" should goes to infinity we have A — K~ 1 e 2 @ N (Fi(f3) + Fi{J3)) where 
Fi([3) — e 2/3 — 1 and F 2 {(3) = (Q — 2) e t$_ x ■ So the action is (apart from a 
positive diverging term) sum of two positive terms, one increasing (with respect 
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to (3) and one decreasing, as it is seen in the figure (JIJ. The minimum value is 
achieved when one puts this terms nearly equal. 

The minimization of the action, therefore is done by letting 4 sinh 2 (3 ~ 
(fi-2) leading to j3 = ±(3 with /? = § cosh" 1 f . So the solution is Q n = e ±ft >™, 
which is again the exact solution 3 . 

As a result of introducing a cut-off, the minimization depends on it. That 

is, if we translate the cut-off point from N to N + c, the minimum point moves 

from the point where the derivative of F\ + F 2 is zero to the point where the 

derivative of F\ + F2 equals to — 2c(F\ + F 2 ). The solution is to renormalize 

the regularized action via introducing the physical gravitational constant as: 
K -i = K -l e 2/3AT 

physical 



7 Case: Q ^ 2 and 5„ ^ 

Now the action would be: 

oo 

A= £ ( K - 1 |g„-Q„_ 2 | 2 + ( K - 1 (r!-2) + 5„-i)|Q„-i| 2 ) (27) 

n— — oc 



Guessing a solution of type e^ n ' 5 we have: 

oo oo 

first term = kT 1 ^ | e ff(«) _ e ff(™-2)|2 ^ K -l ^ 4 sinh 2 .a'e 9+9 * (28) 



n— — oo n— — oo 



second term = ^ (K -1 (fi - 2) - S n )e 9+3 * (29) 

n— — oo 

Again we make the action extremum by letting these two terms equal. To 
see why this is true, we note that the Lagrangian (|15|l consists of two terms 
Ci = k _1 |Q„+i - Qn-i\ 2 and C 2 = k _1 (0 - 2)|Q„| 2 + S n \Q n \ 2 (where we have 
shifted n by one). On the equation of motion i|16|) the second term would be 
C 2 = K _1 <3^(Qn+2+Q„-2)-2K _1 |Q„| 2 . Expanding g (nil) andg(n±2) around 
g(n) we have d ~ K~ 1 e 9+ ^ \e 3 ' - e~ a '| 2 and C 2 ~ K~ l e B+g * {e 2g ' + e~ 2 f' - 2). 
These are equal provided we assume either g is real or purely imaginary. Note 
that this is an approximate result since we expanded g and choose only the first 
term of expansions. 

The condition that these two terms be of the same order leads to: 

g = ± J dnsinh- 1 n ~ 2 ~ kS " (30) 

which is clear that is either real or purely imaginary as we assumed. Since we 
had expanded the first term, we have to have g" « g' for usability of this 
solution. For S„ = n a with a an arbitrary number, this condition holds. 
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8 cosmological constant 

As an example let us find the solution for the special case in which the matter 
contribution is cosmological constant. The source term is: 

An M»l 

AM 1 ^ \n\ ~ 1 1 ' 

where A is proportional to the cosmological constant. 

• |n| ~ 1 In this case we have: 

g(n) = ±^-^(Sl - 2 + K A|n| 1 /2)3/2 (2(r! _ 2 ) _ s K A\ n \^/ 2 ) 

„ 4(n-2)'/» (Si- 2)V» 
" 15A2 K 2 ± 2 M 
which is very close to exact solution [5] . 

• \n\ » 1 In this case we have: 



8 l/rj-2 + KAn\ , _i Sl-2 + nAn 
K n ) = TT S o a smh 



kA 2 V 4 / V 4 



1 /fi-2 + KAn / 0-2 + rtAn 1 , , SI -2 + nAn \ 
-IV 4 V 1 + 4 + 4 Smh V 4 ) 

(33) 

Using the asymptotic forms of Bessel functions one can see that this so- 
lution is in agreement with the exact ones, which are J n /2+i/A(l/A) and 
Y n/2+1/A (l/A)®. 

9 Other matters 

In general, matter sources in large or small n limit lead to a source term like 
S n — An". For different values of a the integral for g(n) can be calculated: 

• a = — 2 (large n, massless matter) 

n 

9{n) = 



((SI 2 -3SI + 2) n 4 + kA(2S1 - 3)n 2 + kA 2 ) 



((SI 2 - 3fi + 2) n 4 + kA(2Q - 3)n 2 + kA 2 ) shah -1 (J^+Sl-2 



, (Sl-2)n 2 + KA (Sl-l)n 2 + nA kA „ kA 
- ! in; \ \ I ± ^— d ± ^— \ / — + fi - 2\ I — + < > - I ■ 
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where F is the elliptic function defined as: F[^>\m) — J Q d6(\— msin 2 #)~ 1//2 . 

• a = 1/2 (small n, massless matter with the quantum ambiguity [i = 0) is 
the same as small n, cosmological constant. 

• a = 1 (small n, massless matter with = 1/4) is the same as large n, 
cosmological constant. 

10 Continuous Approximation 

The continuous approximation can be simply achieved by expanding any Q 
around Q n up to second order. The equation of motion is 

4Q"(n) = («5(n) + 2)Q(n) + 0(3) (35) 

and the conserved quantity emerging from the phase invariance is 

J = 3(Q(n)Q*'(n)) + 0(3) (36) 

This is the conserved current expected from the above equation of motion. 

11 Conclusion 

We saw that it is possible to study the LQC difference equation using action 
functional formulation. We have shown that one can get acceptable solutions 
by guessing one and making the action optimized. The results are in agreement 
with numerical and other analytical methods in the literature (see, e.g. [SI 
IHlEIj)- Therefore the method presented here is useful for obtaining analytical 
solutions in different regimes and for different matter sources. 
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